New High Field State of Flux Line Lattice in Unconventional Superconductor CeCoIns 
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Ultrasound velocity measurements of the unconventional superconductor CeCoIns with extremely 
large Pauli paramagnetic susceptibility reveal an unusual structural transformation of the flux line 
lattice (FLL) in the vicinity of the upper critical field. The transition field coincides with that at 
which heat capacity measurements reveal a second order phase transition. The lowering of the sound 
velocity at the transition is consistent with the collapse of the FLL tilt modulus and a crossover to 
quasi two-dimensional FLL pinning. These results provide a strong evidence that the high field state 
is the Fulde-Ferrel-Larkin-Ovchinikov phase, in which the order parameter is spatially modulated 
and has planar nodes aligned perpendicular to the vortices. 

PACS numbers: 



In almost all superconductors, once the energy gap in 
the electronic spectrum opens at the critical temperature 
Tc, only the gap amplitude, but not the shape and sym- 
metry around the Fermi surface, changes within the SC 
state. The only exceptions to this well-known robust- 
ness of the SC gap symmetry have, until now, been re- 
ported for UPt3[l], Sr2Ru04[2], PrOs4Sbi2[3], as well as 
for superfluid ^He. In these materials, multiple SC phases 
with different symmetries manifest themselves below Tc. 
Moreover, in these superconductors the SC order param- 
eter possesses a multiplicity, with a near-degeneracy of 
order parameters with different symmetries. Tuning the 
pairing interaction by an external perturbation such as a 
magnetic field then causes a SC state of given symmetry 
to undergo a transition to a different SC state. 

On the other hand, in the early 1960's, Fulde and 
Ferrell and Larkin and Ovchinikov (FFLO) have de- 
veloped theories [4], different from the above-mentioned 
mechanism, for a novel SC phase with a different SC 
gap function. Generally, in spin singlet superconduc- 
tors, superconductivity is suppressed by a magnetic field 
as a consequence of its coupling to the conduction elec- 
tron spins (Pauli paramagnetism) or to the orbital an- 
gular momentum (vortices), both of which break up the 
Cooper pairs. The FFLO phase appears when Pauli pair- 
breaking dominates over the orbital effect [5-9]. In the 
FFLO state, pair-breaking arising from the Pauli effect 
is reduced by forming a new pairing state (/cj, —k + ql) 
with \q\ ~ 2iibH /hvp {vf is the Fermi velocity) be- 
tween exchange-split parts of the Fermi surface, instead 
of (/cj, — /c|)-pairing in ordinary superconductors. As 
a result, a new SC state with different order parameter 
appears in the vicinity of the upper critical field i^c2- 
One of the most fascinating aspects of the FFLO state 
is that Cooper pairs have a drift velocity in the direc- 
tion of the applied field, and develop a spatially modu- 
lated order parameter and spin polarization with a wave 
length of the order of 27r/|q|. In spite of enormous ef- 



forts to find the FFLO state, it has never been observed 
in conventional superconductors. The question of ob- 
serving the FFLO phase in an unconventional supercon- 
ductor has only been addressed more recently. In the 
last decade, heavy fermion superconductors CeRu2 and 
UPd2Al3 have been proposed as candidates for the ob- 
servation of the FFLO state, but subsequent research has 
called the interpretation of the data for these materials 
in terms of an FFLO state into question[10]. 

Recently it was reported that CeCoIns is a new type 
of heavy fermion superconductor with quasi two dimen- 
sional (2D) electronic structure and an effective elec- 
tron mass m* ~ 100me[ll]. CeCoIns has the high- 
est transition temperature (Tc=2.3 K) of all heavy 
fermion superconductors discovered until now. Subse- 
quent measurements have identified CeCoIns as an un- 
conventional superconductor with, most likely, d-wave 
gap symmetry [12, 13]. Very recent reports of heat ca- 
pacity measurements in field parallel to the a6-plane of 
CeCoIns have raised great interest, because they pos- 
sibly point to the occurence of a FFLO phase[14, 15]. 
The phase transition from SC to normal metal at the 
upper critical field parallel to a6-plane, i^j2 5 is of first 
order below approximately 1.3 K, in contrast to the sec- 
ond order transition in other superconductors [13-16]. A 
second order phase transition line, which branches from 
the first order H^^2~^^^^ decreases with decreasing 
T, was observed within the SC state below 0.35 K[14]. 
CeCoIns satisfies the requirements for the formation of 
FFLO state. First, it is in the very clean regime, i.e. 
the quasiparticle mean free path is much larger than the 
coherence length ^. Second, the Ginzburg-Landau pa- 
rameter = A/^ 10, with A the penetration length, is 

very large. Thirdly, i^j2(^ — ^) extremely high 

value of 11.9 T, owing to a very large conduction elec- 
tron mass and to the two dimensionality. This situation 
is favorable for the occurence of the FFLO state because 
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the Pauli effect may overcome the orbital effect. Fourth, 
the pairing symmetry is most hkely to be (i-wave, which 
greatly extends the stability of the FFLO state with re- 
spect to a conventional superconductor [8]. On the basis 
of some of these arguments, the possibility of a high-field 
FFLO state has been evoked[14, 15]. However, there still 
is no corroborated evidence for the FFLO state, because 
of other possible origins for a second order transition, 
such as some sort of magnetic transition in the heavy 
fermion system. Therefore, a more detailed experimen- 
tal investigation of the vortex state is required in order 
to shed light on this subject. 

The most salient feature of the FFLO state predicted 
by a number of theories is that the SC order parameter 
has planar nodes aligned perpendicularly to the applied 
field. In particular, as pointed out in Ref. 6, the fiux 
lines are divided into segments with distance ~ 27r/ | q|. 
Therefore, in order to establish the existence of a possi- 
ble FFLO state, it is particularly important to elucidate 
the structure of the fiux line lattice (FLL), which in turn 
is intimately related to the electronic structure. Here, 
we present ultrasound velocity measurements which pro- 
vide direct information on the FLL structure in CeCoIns. 
In ultrasound experiments, sound waves are coupled to 
the FLL when the latter is pinned by crystal lattice de- 
fects. The ensuing modified sound dispersion allows one 
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FIG. 1: The transverse sound velocity Vt as a function of 
temperature in two different configurations, (a) the Lorentz 
mode H II k II [100], u ||[010] and (b) the non-Lorentz-mode, 
u II H II [100], k II [010]. The configuration (a) corresponds to 
a flux line bending mode. The inset to Fig. 1(b) shows in 
zero field and in the normal state above i/'L. 



to extract detailed information about the FLL-crystal 
coupling. We find a distinct anomaly of the sound veloc- 
ity, the nature and magnitude of which provides strong 
evidence for a vortex state with modulated order param- 
eter as predicted for the FFLO phase. 

The ultrasound velocity measurements were performed 
on high quality single crystals of CeCoIus with tetrago- 
nal symmetry, grown by the self-fiux method. Ultrasonic 
waves, with a frequency of 90 MHz, were generated by 
a LiNbOs transducer glued on the surface of the crystal. 
A pulse echo method with a phase comparison technique 
was used for the sound velocity measurements. The rel- 
ative change of the sound velocity was measured by the 
phase change of the detected signal. The resolution of 
the relative velocity measurements was about 1 part in 
10^. In all measurements, the magnetic field H was ap- 
plied parallel to the 2D planes (H || [100]). The inset to 
Fig. 1(b) shows in zero field, and in the normal state 
above i^j2- zero field, shows a steep increase be- 
low Tc, indicating a stiffening of the crystal lattice at the 
superconducting transition. In all experiments presented 
here, the crystal lattice response is limited by its shear 
modulus, Cqq. 

Experiments on the FLL were always carried out in the 
field-cooled condition. Figures 1 (a) and 1(b) display the 
transverse sound velocity as a function of temperature 
in two different configurations with respect to the polar- 
ization vector u, the sound propagation vector k, and H. 
We distinguish between the Lorentz force mode (Fig. 1 
(a), u ± H) where the sound wave couples to the vortices 
through the induced Lorentz force Fl ^ A~^(u x H) x B, 
and the non-Lorentz mode (Fig. 1 (b), u || H) where the 
sound wave couples only to the crystal lattice. As shown 
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FIG. 2: The relative shift of the transverse sound veloc- 
ity Avt /vt as function of H at 100 mK for the Lorentz 
(H ± u) and the non-Lorentz modes (H || u). The inset de- 
picts the difference between Lorentz and non-Lorentz modes, 
Avt = Vt{u ± H) — v°(u II H), which can be attributed to 
the contribution of the FLL. The dashed line shows a fit to 
Eq. (1). 
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in Figs. 1 (a) and (b) and Fig. 2, the transverse sound 
velocity in the Lorentz mode, in which flux motion is par- 
allel to the a6-plane, is strongly enhanced compared to 
that in the non-Lorentz mode. This indicates that the 
transverse ultrasound strongly couples to the FLL in the 
Lorentz mode. The difference between the sound veloc- 
ity in the Lorentz mode and in the non-Lorentz mode, 
Avt = v^{u 1. H) — v^{u II H), can be regarded as being 
the contribution of the FLL to the sound velocity. In 
the Lorentz mode. Fig. 1(a), the transverse polarization 
of the sound has the FLL undergo a long-wavelength ( 
I k I A <C 1 ) bending mode, limited by the FLL tilt mod- 
ulus in the limit | k |-^ 0, c{4(| k |^ 0) = ^Vmc The FLL 
contribution can then be written as [17] 



Aiopt^tl + (*oS)i/2|k|VMoic' 

with p the mass density of the superconducting material, 
jc its critical current density, and = h/2e. 

Figure 2 shows the relative shift of the transverse sound 
velocity Av^/v^ as a function of field at very low tem- 
perature. The hysteresis arising from the trapped field 
was very small. Below approximately 9 T, Av^ /v^ for the 
non-Lorentz mode is nearly i^-independent while Av^ /v^ 
increases steeply with H for the Lorentz mode. By fitting 
the increase to Eq. (1), we obtain a (field-independent) 
critical current density, jc = 7.0 x lO^Am"^. In both 
modes, Av^/v^ decreases gradually at higher field with 
i7, followed by a jumplike decrease at i^j2- This sharp 
drop, which occurs in a very narrow i^-range AH/H^^2 < 
1%, is an indication of the first order transition, and 
disappears at high temperature where the transition be- 
comes second order. 

In Figure 3, several curves of the relative shift of the 
FLL contribution Avt / Vt below Hc2 are shown as a func- 
tion of temperature. As the temperature is lowered be- 
low the transition temperature in magnetic field Tc(i^), 
Avt/vt starts to increase. At fields of 9.5 T and higher, 
Avt/vt changes its slope with a distinct cusp at T*, as 
marked by the arrows, while at H=8.5 T, no anomaly is 
observed. Obviously, T* separates two different vortex 
states with different transverse sound velocities; Avt/vt 
below T* is smaller than the Avt/vt extrapolated from 
above T* . This indicates a sudden decrease of FLL cou- 
pling to the crystal lattice below T* . Such an anomaly in 
the ultrasound velocity was not observed when the mag- 
netic field was swept, probably because it occurs near the 
maximum of Avt{H). 

Figure 4 displays the H—T phase diagram below 0.7 K. 
In this temperature range, the transition at Hc2 is of first 
order as indicated by the sharp drop of the ultrasound 
velocity (see Fig. 2). In the same figure, we plot the Hc2- 
line reported in Ref. 14. Although H^^2 crystal is 

slightly higher than reported in Ref. 14, T* seems to 
coincide well with the reported second order transition 
line[14, 15]. On the basis of our results, we conclude that 
the second order phase transition is characterized by a 
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FIG. 3: The relative shift of the transverse sound velocity 
arising from the contribution of the FLL, Avt/vt, as a function 
of temperature, below the upper critical field ifj2- Avt/vt is 
obtained as the difference between the non-Lorentz mode and 
the Lorentz mode sound velocities. For clarity, the zero levels 
of the different curves are vertically shifted. At T*, marked by 
arrows, Avt/vt exhibits a distinct cusp. The dashed lines are 
Avt/vt extrapolated from above T*. Dotted arrows indicate 
Tc{H). 



sudden change of flux line pinning, due to a structural 
transition of the FLL. The present results definitely rule 
out the possibility that some kind of magnetic transition, 
or a structural phase transition of the crystal lattice, is at 
the origin of the transition. This can be seen by looking 
at Av^ /v^ in the non-Lorentz modes (Fig. 1 (b)), which 
show no anomaly at T*. 

We here discuss several possible origins for a struc- 
tural change of the FLL. A reduction of the ultrasound 
velocity with H at the transition would occur when the 
vortex lattice melts into a liquid, as observed in layered 
high- Tc cuprates and amorphous thin films [18]. However 
we can completely rule out this scenario, because the en- 
tropy jump at T* reported in Ref. 14, 15 is significantly 
larger than that expected for vortex lattice melting. In 
addition, flux pinning vanishes at a FLL melting transi- 
tion, causing Avt/vt to abruptly drop to zero, at odds 
with the experimental observation. A second possibility 
is a FLL order-disorder transition in the presence of pin- 
ning. However, this should be accompanied by an abrupt 
increase of the measured critical current and of the ultra- 
sound velocity, whereas the present data show a decrease. 
Such an effect has been discussed in the light of a possible 
FFLO state in CeRu2 and UPd2Al3[19]. Thirdly, there 
might be a possibility that the nature of FLL pinning 
changes above the transition. For instance, new pinning 
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FIG. 4: Experimental H — T phase diagram for CeCoIns be- 
low 0.7 K for (H II [100]). The transition temperature T*, 
indicated by red squares, were obtained from ultrasonic mea- 
surements as indicated by the arrows in Fig. 3. Blue circles 
and solid line depict the upper critical field Hc2 determined 
by the ultrasonic experiments. In this temperature regime, 
the transition at H^^2 is of first order. The dash- dotted line 
is the second order transition line determined by the heat ca- 
pacity reported in Ref. [14]. The region shown by green is in 
the FFLO state. The schematic figures are sketches of a flux 
line above and below the transition. 



centers, such as localized spins, could be induced by the 
transition. This effect may cause the change of the sound 
velocity by stiffening the flux lines. However, such an ex- 
trinsic origin is quite unlikely because of the absence of an 
anomaly in the non-Lorentz mode, as discussed before. 
Finally, it may be that the FLL undergoes a symmetry 
change, e.g. from hexagonal to square, with a concomit- 
tant lowering of the vortex lattice shear modulus Cgg and 
increase of the pinning force. Such a transition has in- 
deed been observed in CeCoIn5[20], but at /j^qH = 0.6 T, 
much lower than the fields involved here. 

The failure of all the above scenarios leads us to con- 
sider the possibility that the anomaly in the FLL pin- 
ning force and the sound velocity is intimately related 
to a change of the quasiparticle structure in the high 
field SC phase, with a reduced electron correlation along 
the direction of the FLL in the high field phase. Start- 
ing from the hypothesis of a transition from a "contin- 
uous" superconductor to one with an order parameter 
structure, modulated in the field direction, with no elec- 
tronic correlations between layers, we can understand 
the modest decrease of Avt/vt. From the critical cur- 



rent density above T*, we deduce the pinning strength 
W of crystalline defects such as small dislocation loops or 
stacking faults. For a three-dimensional superconductor, 
this should be evaluated in the single vortex limit [21]: 

W = {^l^^B^/^fJ^TT^fioXY^^ = 5 X lO-^N^m^ At 
the transition to the phase with modulated order param- 
eter, the vortex tilt stiffness in the limit of the small 
scale distorsions caused by pinning collapses. Then, 
the critical current is given by the expression for a 
2D layer of thickness d, jl^ = W/^^^ B^^'^cl^d. At 
the transition, j^^ must be equal to the experimental 
value 7 x lO^Am"^, an equality, that, given the exper- 
imental value of W ^ is satisfied for a layer thickness of 
4 X 10~^ m. This is remarkably close to the modulation 
27r/ |q|= {me/m*){hkF/eB) = 3.5 x 10"^ m expected for 
the order parameter structure associated with the forma- 
tion of the FFLO state proposed by Ref. 5-9 (the Fermi 
wavevector ^ S x 10^m~^). In the FFLO state, the 
order parameter performs one-dimensional spatial mod- 
ulations with a wavelength of the order of 27r/ |q|, which 
is comparable to ^, along the magnetic field, forming pla- 
nar nodes that are periodically aligned perpendicularly 
to the flux lines. The occurrence of planar nodes leads 
to a segmentation of the flux lines. A schematic figure of 
this state is shown in Fig. 4. This field-induced layered 
structure resembles the vortex state of layered cuprates 
and organic superconductors in magnetic field perpen- 
dicular to the conducting planes. The segmentation of 
vortex lines to a length 27r/ |q| much smaller than the 
characteristic "Larkin length" Lc = {cl^B /^qW^^^)^^^ 
on which they are pinned, means that the vortex part of 
the ultrasound velocity is now governed by the pinning 
of the quasi 2D vortex layers of thickness 27r/| q|. Note 
that a decrease of Avt/vt at the transition can only be 
obtained by a true dimensionality change of the FLL at 
the transition, and not by a mere softening of the FLL 
elastic constants. 

In summary, ultrasound velocity measurements of 
CeCoIns in the vicinity of the upper critical field Hc2 
reveal an unusual structural transformation of the FLL 
at the second order phase transition within the SC state. 
This transformation is most likely characterized by the 
collapse of the FLL tilt modulus and a transition to a 
quasi-2D state. These results are a strong indication that 
the high field superconducting state is the FFLO phase, 
in which the order parameter is spatially modulated and 
has planar nodes aligned perpendicularly to the FLL. 
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